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ABSTRACT 

We present a sample of II M3I Cepheids in stellar clusters, derived from the overlap of the 
Panchromatic Hubble Andromeda Treasury (PHAT) cluster catalog and the Pan-STARRSl (PSl) 
disk Cepheid catalog. After identifying the PSl Cepheids in the HST catalog, we calibrate the PSi 
mean magnitudes using the higher resolution HST photometry, revealing up to 1 magnitude offsets due 
to crowding effects in the ground-based catalog. We measure ages of the clusters by performing single 
stellar population fits to their color-magnitude diagrams (CMDs) excluding their Cepheids. From 
these cluster age measurements, we derive an empirical period-age relation which agrees well with the 
existing literature values. By confirming this relation for M31 Cepheids, we justify its application in 
high-precision pointwise age estimation across M31. 

Subject headings: Cepheids — galaxies: individual (M31) — galaxies: star clusters: general 


1. INTRODUCTION 

Cepheids are young pulsating variable stars, widely- 
known for their utility in distance measurement. They 
exhibit a tight period-luminosity (P L) relation (the Leav¬ 
itt Law, Leavitt fc Pickering|1912 1, a product of the fact 
that this form ot pulsation only occurs for stars in an 
unstable evolutionary phase occupying a thin strip in 
luminosity-effective temperature space. Cepheids also 
exhibit a period-age relation, as more massive stars 
will evolve through the Cepheid strip earlier and at 
brighter magnitudes (Kippenhahn & Smith|1969|. Thus, 
Cepheids can act as standard candies as well as useful 
probes of star formation history in the local universe. 

Testing the period-age relation directly requires an 
age estimate independent of the Cepheids. These ages 
are generally obtained through color-magnitude diagram 
(CMD) fitting applied to a single-age stellar population 
(SSP), which limits rigorous testing to well-resolved stars 
with known distances that are close enough for the main 
sequence turnoff to be visible. Thus, previous tests of 
period-age relations have focused on Cepheid samples in 
our Galaxy or in the Small and Large Magellanic Clouds 
(MCs: SMC and LMC, respecti vely; see for instance, 
Bono et al. 2005 Efremov 20031. The MCs are well- 
resolved in ground-based surveys and have large samples 
of Cepheids and stellar clusters, making them well-suited 
for CMD age estimation. However, they have signifi¬ 
cantly lower metallicity than the more massive galaxies 
in the Local Group where Cepheid distance and age mea¬ 
surements are often applied. 
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The Andromeda galaxy (M31) provides an excellent 
opportunity to constrain the period-age relation at high 
metallicity. This spiral galaxy is close enough that indi¬ 
vidual bright stars in it are well-resolved by the Hubble 
Space Telescope (HST), even within stellar cluster envi¬ 
ronments. Existing period-age relations have bee n used 
to infer disk sta r formation in M31 (for instance, Mag- 
nier et al.||1997|), but not yet directly calibrated there. 


Thanks to dedicated CCD surveys, large samples of M31 
Cepheids have been obtained i n recent years . The DI¬ 
RECT project, summarized by Macri (20041, searched 
for Cepheids and de tached eclipsing binari es in a num¬ 
ber of M31 fields. IVilardell et al. (2007) identified a 


sample of 416 Cepheids m the process of searching for 
eclipsing binaries, with photometry sufficient to charac¬ 
terize the pulsation modes of 356 Cepheids. The largest 


M31 Cephei d sample to-date was obtained by [Kodric 


et al. (2013), using the Pan-S TARRS 1 (PSl) survey 
telescope PAndromeda dataset (Lee et al. 2012). This 
sample avoids many of the issues which plagued earlier 
searches, namely: limited spatial coverage, few epochs, 
and short observational campaigns (which can introduce 
bias against finding Cepheids in certain period ranges). 
The PSl Cepheids were found via single-camera differ¬ 
ential photometry covering the entire M31 disk. The 
photometry were obtained in 183 epochs, spanning a 
year. The lightcurves are well-sampled, and allow for 
robust period determination and Cepheid characteriza¬ 
tion. A total of 1440 of these Cepheids were identified as 
fundamental-mode pulsators. 

Progress has also been made on individual-star pho¬ 
tometry in M31. The Panchromatic Hubble Andromeda 


Treasury, o r PHAT (Dalcanton et al. 2012 Williams 


et al. 2014), is a large survey ot a third of the M31 


disk in 6 hlters with HST. This survey provides another 
means to study individual Cepheids found using high- 
caden c e ground-based photometry . For instance, IRiessI 
et al.| (20121), Kodric et al. (2015[), and Wagner-Kaiseij 


et al. 


the NIR period-luminosity relationship and improve the 
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determination of the M31 distance modulus. PHAT ob¬ 
servations have also enabled the develop ment of a large 
and well-characterized star cluster sample (Johnson et al. 


20151, most of which are well-resolved in the HS'i' data. 

In this paper, we measure ages for individual Cepheids 
by identifying them in stellar clusters with secure mea- 
sured ages. Specifically, we cross-correlate the IKodric 


et al. 


(2013) Cepheid catalog with the Johnson et al. 
star cluster catalog. The final sanipfe of 11 high- 
quality candidates adds to the approxima.tely 2 3 Galactic 
open cluster Cepheids (|Anderson et al.||2013D, together 
providing a solar metaln city complement to samples in 
the LMC and SMC (e.g., Pietrzynski & Udalski 1999). 

In Section we describe the PHAi' and Psl data 
products which we utilize. Section describes the iden¬ 
tification of the Cepheids in the PHAT data, the process 
and results of dephasing the PSl magnitudes and the de¬ 
termination of isochrone fits for the clusters. In Section 
l^we compare these cluster age determinations with sev¬ 
eral period-age relations in the literature. We conclude 
with a brief summary in Section 

2. OBSERVATIONAL DATA 

2.1. PAndromeda Cepheids 

We select Cephei ds from the Pan-ST ARRS Cepheid 
sample presented by Kodric et al. (2013), which is large 
and well-characterizeol Observations were made in the 
ipi and rpi filters over 183 epochs from 2010 July 23 
to 2011 August 12 (with a half-year gap in the sched¬ 
ule). Before inclusion in the final sample, the candi¬ 
dates passed cuts base d on CMP loca tion and lightcurve 
Fourier parameters. [Kodric et al.| manually classi¬ 
fied a subset of the lightcurves and defined regions in 
Fourier space corresponding to fundamental mode, first- 
overtone, and Type-II Cepheids based on these classifi¬ 
cations. This enabled them to automatically classify all 
of the objects in their catalog. 

The PSl lightcurves are well-sampled, allowing us to 
extend them to the epoch of HST observation with min¬ 
imal error. We can thus correct for the phase of the 
Cepheid at the times of our HST observations, and derive 
accurate mean magnitudes from even the small number 
of HST observations. The process by which we extend 
the PSl lightcurves is described in Section 

2.2. PHAT Clusters 

The most reliable way to assign ages to Cepheids is 
to associate them with a single age stellar population. 
At the distance of M31, the main sequence turnoff for 
even young populations is near the HST detection limit. 
Though it is possible to assign ages to recent burs ts of 
star formation in any region of the M31 disk (e.g., Jen¬ 
nings et al. 2012), the results suffer from potential bias 
due to the presence of overlapping populations and small 
numbers of coeval stars. We are interested in robustly 
testing the period-age relationship, and bias in the age 
estimates (especially from older underlying populations) 
could signihcantly affect the result. Thus, we focus on 
Cepheids near high-conhdence clusters, where we can ex¬ 
pect to obtain strong age constraints. 

The Andromeda Project (hereafter AP) provides an 
M31 cluster sample of unprecede nted size and quality , 
derived from the PHAT survey (Johnson et al. |2015). 


The PHAT data consists of HST Wide Field Camera 
3 (WFC3) and Advanced Camera for Surveys (ACS) 
imaging in 6 hl ters covering ^ 0.5 de g^ of the north¬ 
east M31 disk (jDalcanton et al. 2012). The AP clus¬ 
ters were identified through visual inspection of optical 
ACS/WFC F475W and F814W search images by citizen 
scientist volunteers. We use the hnal AP cluster catalog 
of 2,753 clusters for this work. 

2.3. Cross-Matched Cluster Cepheid Sample 

To form our sample of candidate cluster Cepheids, we 
first search for PAndromeda Cepheids whose reported 
positions place them within the extent of an AP clus¬ 
ter aperture. Since the AP cluster radii and centers are 
somewhat uncertain, we first define the extent generously 
by expanding the AP radius by a factor of 1.2; however, 
all matches found were within the nominal cluster radii. 
This search yields 9 fundamental mode (FM) Cepheid s 
and 2 left unidentified (UN) by Kodric et al. (2013). 
Fundamental, first-overtone, and I'ype-H Cepheids must 
be studied independently for precise period-age-relations; 
but since all but two of our objects are firmly identified 
as FM pulsators and our sample size is already small, 
we proceed under the assumption that they are all FM 
Cepheids. Table details the basic properties of this 
set. The spatial and period-luminosity distribution of 
the sample is depicted in Figures and In addi¬ 
tion, we note that all of the clusters sele ctea have final 
weigh ted cluster fractions (discussed in [Johnson et al. 
2015) fcist,w > 0-6, suggesting that all of these Cepheids 
tall within highly-probable clusters. 

The PHAT survey implementation can be exploited to 
extract multi-epoch data for most points in the survey 
grid. This is due to the fact that the HST tiling is based 
on WFC3/IR, which has a substantially smaller held of 
view than ACS/WFC; so simultaneous pointings of these 
cameras generate large overlap in the ACS images. This 
means that for all but one of the Cepheids we have deep 
HST imagery from at least two different times, and thus 
a means to assess variability from the HST data alone. 
This extra information can be used to distinguish the 
Cepheid from its neighbors at very high conhdence. The 
number of HST images available for the clusters under 
consideration ranges from 1 to 4 with a mean of 2.7. 

3. ANALYSIS 

Before we can study them in detail, we must conh- 
dently id enti fy the Cepheids in the crowded PHAT helds 
(Section 3.1). This enables us to improve the Cepheid 
mean m agni tudes using the superior resolution of HST 
(Section 3.2); and to d erive independent age estimates for 
the Cepheids (Section 3.3). The PSl positions are par¬ 
ticularly uncertain in clusters which HST resolves into 
multiple bright stars. In order to select the Cepheids in 
the HST data robustly, we examine both inst ability strip 
posit ion (using theoretical boundaries from Bono et al.j 
2005) and PHAT field-to-held variability. 


3.1. Cepheid Identification in the HST Data 

We hrst attempt to identify each Cepheid as a resolved 
source in the HST photometry. Spatial agreement is 
checked by cluster proximity. As the cluster helds are 
crowded and the PS-1 determined positions are given to 
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TABLE 1 

PHAT/PS 1 Cluster Cepheids 


CC ID 

AP ID 

PS ID 

ClusterFrac 

RA (J2000) 

DEC (J2000) 

mode 

Period, rp\ (days) 

ipi,HST — ipi.PSl 

1 

477 

PSO-JOll.4286-1-41.9275 

0.9881 

11.42865 

41.92759 

EM 

4.582 

1.19 ±0.08 

2 

1539 

PSO-JOll.4279-1-41.8642 

0.9358 

11.42797 

41.86422 

UN 

6.08 

0.99 ±0.21 

3 

3928 

PSO-JOll.6227-1-41.9637 

0.8738 

11.62273 

41.96373 

FM 

6.213 

0.09 ±0.06 

4 

2831 

PSO-JOll.6374-1-42.1393 

0.9059 

11.63742 

42.13933 

UN 

7.928 

0.88 ±0.05 

5 

3050 

PSO-JOlO.9769-1-41.6171 

0.8461 

10.97692 

41.61718 

FM 

8.829 

0.20 ±0.10 

6 

5216 

PSO-J011.0696-I-41.8571 

0.7049 

11.0696 

41.85714 

FM 

14.353 

0.48 ±0.22 

7 

2113 

PSO-JOll.6519-1-42.1286 

0.9807 

11.65193 

42.12865 

FM 

15.429 

0.11 ± 0.13 

8 

2587 

PSO-JOll.7119-1-42.0449 

0.8284 

11.71191 

42.04491 

FM 

17.18 

0.29 ±0.10 

9 

2967 

PSO-J011.0209-I-41.3162 

0.6609 

11.02092 

41.31625 

FM 

19.566 

0.65 ±0.10 

10 

1082 

PSO-JOll.2797-1-41.6217 

0.9087 

11.27976 

41.62174 

FM 

26.499 

0.31 ± 0.03 

11 

1540 

PSO-JOll.3077-1-41.8500 

0.984 

11.30773 

41.85003 

FM 

35.75 

0.24 ± 0.06 


Note. — Tablelists the C luster Cepheids (CCs) identified from the PSl/PHAT data, and provides some basic information about each. The AP 
ID identifies the cluster in the [Johnson et al.|||2015[ catalog and fcist,w represents the final cluster quality parameter discussed there. The PS ID 
identifies the Cepheid from th e Pb catalog, i'he pe riod in the rpi filter; Cepheid RA and Dec; the mode classification (FM for fundamental mode, UN 
for unidentifed) determined by|Kodric et al.|(|2013|l are also provided. Finally, the last two columns indicate the HST-determined correction to the PSl 
mean magnitudes, as discussed in Section 



Fig. 1.— M31 overview; the full PHAT ACS coverage is outlined 
in purple, AP cluster candidates are purple diamonds, and the 
Kodric Cepheids are marked by red squares. The cluster Cepheids 
we located are denoted by large blue circles. Background image 
credit: Robert Gendler 


an accuracy of only 0.36", we consider all stars within 
the AP cluster radius as candidates. We consider two 
additional pieces of evidence when cross matching the 
Cepheids; location on the CMD and variability. 

Proximity to the instability strip: We adopt the “non- 
canonical”, solar-metallic it y inst ability strip boundaries 
derived by Bono et al. (2005) for fundamental pul- 
sators satisfying 0.17 < log^Q(period/days) < 2.02. 
In order to transform the boundaries from theoreti¬ 
cal TpfF-Lumi nosity space to HST filter space, we use 
the MATCH (Dolphin 2002) fake utility to sa mple the 
Padova isochrones (discussed in Section 3.3) at vari¬ 
ous ages and thereby populate the instability strip with 
model stars. The model stars which fall between the blue 
and red edges and within the valid period range (deter- 
mined using the pulsational relations derived in |Bono| 


et al. 2000) are then plotted in HST F475W-F814W7 
F814W color-magnitude space. The quadrilateral which 
most-tightly encloses these points in CMD space defines 
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24 



10 ^ 

Period (days) 


10 ^ 


Fig. 2.— The Kodric fundamental-mode Cepheid sample, with 
the cluster Cepheid candidates in black, in period-Wesenheit space. 
The Wesenheit index is a color-corrected magnitude measurement 
intended to reduce scatter due to differentia l extinction (in thi s 
case, derived using the rpi and ipi bands in [Kodric et al ||2013||. 
We observe that our candidates are fairly representative ot tlie 
Kodric sample, though there is some indication of a bias towards 
brighter indices (as expected from dense cluster crowding). Note 
that while the uncertainties in the pulsational periods are too small 
to be visible on this plot for most of our sample, there are two 
objects with significant values. 


our observational instabili ty strip. Indeed we ha ve used 
the trans formations fro m Sirianni et al. (2005|, as de¬ 
scribed in Girardi et al. ( 20U8|) (as these transformations 
are used by iVIAi'GH to yield filter magnitudes). 

Membership in the final filter-transformed instability 
strip is uncertain due to errors in the reddening adopted 
for each cluster. We estimate Ay by fit ting the CMD of 
the full cluster (as described in Section 13.3 ). The uncer¬ 
tainty in the reddening estimates from these fits is on the 
order of 0.1 mag. To ensure that we catch all reasonable 
candidates, we expand the edges of the strip to cover an 
Ay range of ±1.5 mag around the best fit, as illustrated 
in Figure This is a very large overestimate intended 
to prevent accumulated error from the HST photome¬ 
try, the filter corrections, and pulsational models from 
excluding too many objects from consideration. 

Variability: For each cluster, we identify all ACS ob¬ 
servations from the PHAT survey whose borders contain 










































































4 


Senchyna et al. 


the cluster. We select a reference field in which the clus¬ 
ter is fully contained and that is not too close to the 
chip gap (where distortion potential is highest). We start 
with the St photometry catalogs (described by |Williams 
et al.|2014|) , and apply an additional cut in S /N (greater 
than 4 in both filters). These catalogs are produced by 
the THAT reduction pipeline, and do not include quality 
cuts on PSF shape or crowding; we wish to maximize 
catalog completeness at the expense of photometry qual¬ 
ity. In order to match the stars in different overl apping 
observ ations, we use pattern-matching code from |Groth| 
(1986) to determine RA/Dec offsets from the brightest 
stars in each field. For each catalogued source in the 
reference field, we search a circle of radius .072 arcsec- 
onds in the other images for matches. This radius is on 
the order of the FWHM for ACS, and our results indi¬ 
cate that it matches all sources of interest. If more than 
one source is found in any of these circles, we consider it 
a non-detection in the corresponding field. The sources 
for which we obtain at least 2 exposures are retained for 
analysis. We compare the resulting catalog to the full st 
reference field catalog in the instability strip region to 
ensure that no candidates are missed due to alignment 
issues. 

Before we can compare the PHAT and PSl data, we 
must extend the PSl lightcurves to the HST obse rvation 
epoch s . For each Cepheid in their final sample, [Kodric 


et al. (2013) provide the raw unfolded lightcurves and 
the coefficients of a fifth-order Fourier series fit to the 
folded lightcurves. Unfortunately, they do not appear to 
provide a zero-epoch for these determinations. We use 
a simple scalar minimization routine to determine the 
zero epoch which leads to the best agreement between 
the Fourier series and the data. The result is checked 
by-eye, and the uncertainty in this epoch determination 
is computed by bootstrap resampling. This procedure al¬ 
lows us to compute the Cepheid magnitude expected by 
PSl at an arbitrary epoch. Uncertainties in the PSl ex¬ 
pected magnitudes are propagated by Monte Carlo meth¬ 
ods, assuming that the parameters and associated errors 
provided by Kodric represent Caussian distributions. 

We compute several variability indices inspired by 


Welch & Stetson (1993) for each matched source. These 
quantify the degree to which the magnitudes in two fil¬ 
ters are correlated from field to field. They are calculated 
for our pair of HST ACS magnitudes F475W and F814W 
for each source, as well as for the source F814W and the 
model expected ipi Cepheid magnitude. We denote fil¬ 
ter 1, filter 2 magnitude in the field by Flk,F2k (re¬ 
spectively), the corresponding pipeline-determined Pois¬ 
son uncertainty in each by CFfi,k,<Jf 2 ,k, and the variance- 
weighted mean for each filter by FI. We then define the 
variability index I to be 


I = 


1 


E 


/Fife-FI 


'J'nin - 1) “ V 0-/1. 


F2fc - F2 
crf2,k 


where n is the number of fields in which the source was 
matched. We denote the index computed with the 2 
ACS filters as Iacs; and with the F814W and expected 
ipi as /acs-psi- Note that this index is the covariance of 
the measurements normalized by the photometric uncer¬ 
tainty. For a truly non-variable source with good pho¬ 


tometry and good background subtraction, the variation 
in measured brightness in two filters should be indepen¬ 
dent. Thus the expectation value of I for non-variable 
stars with ideal photometry is zero. A higher value of I 
indicates a higher degree of correlation between bright¬ 
ness changes in the two filters. A high value of the index 
may indicate a truly variable source; but it can also re¬ 
flect issues with photometry, background subtraction, or 
source matching. Here, our situation is simplified - we 
expect that there is a variable within the field being ex¬ 
amined, and the indices are used to help us identify which 
candidate resolved source it is most likely to be. 

We select as candidates the sources whose weighted- 
average ACS magnitudes place them inside of the gener¬ 
ous instability strip for the cluster (yellow-outlined re¬ 
gion in Figure |^. We identify the putative Cepheid 
as the candidate source with the maximum variability 
index, considering both filters equally. Since the pho¬ 
tometric errors of the ACS observations are substan¬ 
tially smaller than those assigned to the PSl magnitudes 
by our Monte Carlo uncertainty estimation, the ACS 
variability index usually dominates (especially when the 
brightness changes are slight). This selection process is 
illustrated for each cluster in Figure For all clusters 
except AP5216, we find a clear vari^ility standout in 
the expected region of the CMD. Since AP5216 was ob¬ 
served in only one PHAT field, no variability information 
is available for the HST sources; fortunately, the gener¬ 
ous instability strip encloses exactly one cluster source, 
which we identify as CC6. 


3.2. Cepheid Dephasing 

We now proceed to compare the HST and PSl pho¬ 
tometry of the Cepheids. The cluster environments are 
crowded, so we expect the PSl measurements to be sys¬ 
tematically brighter due to source confusion. With our 
Cepheids identified in the HST data, we can test the ac¬ 
tual size of this effect. In this section, we use the absolute 
flux estimate from the HST photometry to re-estimate 
the Cepheid mean magnitudes. This requires that we 
compare the HST photometry to the PSl models in a 
common filter system. 

Since we are only concerned with Cepheids, whose 
SEDs are fairly homogeneous and not unusual, a stan¬ 
dard color-based transformation is appropriate here. 
Both the F814W ACS filter and the ipi filter are close to 
the Landolt (Johnsons/Cousins) I. Tonry et al. (2012) 
has conducted an analysis of the Pan-S i ARKS 1 pho- 
tometric system, which includes transformations to the 
Johnsons/Cousins system derived from a library of 783 
photom etric-standard SEDs. Similarly, Sirianni et al. 
(20051 presents transformations from ACS magnitudes 
to Johnsons/Cousins based on synthetic magnitudes and 
checked with empirical data. We opt to use the quadratic 
relations in both cases. After de-reddening the HST 
data using our cluster Ay, we transform each pair of 
F475W/F814W magnitudes to V/I by iteratively im¬ 
proving the V — I color term. The PSl transformations 
require we use B — V color. Since we are missing a coin¬ 
cident B observation, we opt to estimate the Cepheid’s 
average B — V color using an emp irical period-color rela¬ 
tion from Tammann et al. (2003) for galactic Cepheids. 
To accouiit tor the variation of B — V over the period 
of the Cepheid, we add 0.2 mags in-quadrature to the 
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Fig. 3.— Details of the process of Cepheid identification in the HST data for CC2. Top row: (left) the Kodric lig htcu rve and model 
fit (blue circles and green line, respectively) with filter-transformed and dephased HST data (red triangles; see Section |3.2[ l; (right) same, 
plotted ag ains t time (MJD). Note the large offset between the PSl and HST magnitudes — this difference is discussed in terms of blending 
in Section |3.2| Bottom row: (left) the HST cluster color-magnitude diagram with each star ma rked by a cyan circle. The size of the cyan 
circles scales with the value of the HST F475W-F814W correlation metric discussed in Section |3.1| (stars with very small values are given 
a fixed point size) - candidate instability strip members are additionally marked by red dots, and the putative Cepheid by a red star - the 
instability strip (original, yellow dashed line; and broadened, yellow solid line) are outlined; (right) same as bottom-left, with the HST-fpi 
correlation metric (Section |3.1[ l. Here, we see that one object within the instability strip stands out in both the HST variability metric and 
the HST-ipi correlation metric (the largest cyan circle in both of the bottom plots). The remaining images in this series are presented in 
the Appendix. 


color uncertainty. The reddening is re-appli ed to the fi¬ 
nal ma gnitudes using the prescriptions from |Tonry et ah 
(2012). So for each field in which the Cepheid candi¬ 
date appeared, we obtain an HST-derived estimate of the 
Cepheid brightness in the tpi filter at that field epoch, 
which we can compare to the PSl model magnitude. 

Now we use the HST data to cross-calibrate the PSl 
lightcurves. The PSl magnitude variation provides a 
good estimate of the true change in flux received from 
the Cepheid, but the absolute level of the flux is better 
determined by HST. At each HST epoch, we compute 
the difference in zeropoint-scaled flux between the PSl- 
determined mean magnitude and the PSl model magni¬ 
tude at that time: 


_ 1 n—OAipi 


A/= 10 


- 10 


—0.4ipi „odoi 


then add this to the HST-determined magnitude in the 
ipi filter: 

^PI.HST = -2.51ogio(10-2*^^'"«^/'' + A/). 



Fig. 4.— The difference in 2 pi magnitude between the HST de¬ 
phased means and those determined by PSl, plotted against a 
blending metric. The metric J^CepheidlJ^totai is the ratio of the 
putative Cepheid flux to that of all HST-resolved and cataloged 
sources within 0.86”. The results from each PHAT held are av¬ 
eraged. There is a clear correlation between field crowding and 
brighter PSl magnitudes. 


We repeat this for each field, and propagate the un¬ 
certainties via Monte Carlo methods (assuming Gaus¬ 
sian distributions for the parameters). Thus, for each 
Cepheid, we obtain at least 2 estimates of a calibrated 
mean magnitude in ipi. 

The results of this dephasing step are shown in Fig¬ 
ure]^ We plot the difference between the HST-calibrated 
mean magnitude and the PSl mean magnitude (averag¬ 
ing the results for all PHAT fields the star was located 
in), against a blending measure. The measure is com¬ 
puted by summing the F814W flux of all .dst catalog 


sources within a radius of 0.86” (3.3 parsecs in the disk) 
of the putativ e Cepheid. In de scribing the PAndromeda 
data pipeline, |Lee et al.| (|2012|) reports for one represen¬ 
tative skycell a median FWHM in the 30 best images 
of 0.861". Thus, the metric represents a rough estimate 
of the degree to which blending affects the PSl data. 
The HST-calibrated means are systematically dimmer 
than the PSl means, and there is a trend towards better 
agreement with less excess (HST-cataloged) flux near the 
Cepheid. 

As expected in crowded cluster environments, the PSl 
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16.5 



of convective overshoot. 


1.0 1.2 1.4 

logio(period/days) 

Fig. 5.— The F160W PL relation for our cluster Cepheids, com¬ 
pared with t he empirical linear PL relation found by |Wagner-Kais^ 
|et al.| (|2015|l in M31. The two agree within the uncertainty, which 
IS represented by the filled section; the magnitude errors in our 
data are photometry-derived only (no consideration was given to 
Cepheid phase). Three Cepheids were excluded due to NIR cover¬ 
age issues. Note that the uncertainties in the Cepheid magnitudes 
and periods are smaller than the point size on this plot except for 
two of the period measurements. 

magnitudes appear to be systematically brightened by 
blending. This effect has been studied in-depth pre- 
viously, e.g. in relatio n to th e DIRECT pro ject b; 


Mochejska et al. (2000) and by Vilardell et al. ( 2007| 
In hnd ing a TL relation for their sample, |Kodric et A. 


(|2013|) performs iterative 3 ct clipping to eliminate blends 
and other outliers. We repeat their procedure, with both 
the full set of Kodric FM and UN Cepheids; and with 
the FM Cepheids alone. In both cases, only one Cepheid 
from our cluster Cepheids sample (CCl) is cut (after 5 
iterations in the former, and on the first iteration in the 
latter case). This leaves 2 (4, including the UN) Cepheids 
in the PL-determining set whose PSl measured magni¬ 
tudes are offset at the 0.5-1.0 mag level. 

We also find the Cepheid magnitudes in the near- 
infrared filter F160W. The effects of random-phase ob¬ 
servations are dimini shed substan t ially i n this filter, as 


exploi ted in M31 by iRiess et al. (120^ 


aiser et al. 


Kodric et al. 
i'he bdOOW 


(2015), and Wagner- 
magnitudes ot 8 of our cluster Cepheids (ff were not im¬ 
aged in NIR filters) are plotted in Figure ]^ The empiri¬ 
cal (li near) F160W PL relation derived byfWagner-Kaiser 


let ai.| is also plotted. Our Cepheids show good agreement 
with this relationship within the measured uncertainty. 

3.3. Cluster Age Determination 

We estimate the age of the host stellar clusters by fit¬ 
ting populated isochrones to their CMDs. In particular, 
we use Padova isochrones and assume a Kroupa IMF. 
These isochrones combine the tracks from IGirardi et al.l 
(2000) through the first thermal pulse on the AGB with 
tracks that trea t the TP-AGB in be tter physical detail, 
as described in Mari go et al. (2008). Gorrections from 


Girardi et al. (2010) concerning the AGB are also in- 


corporated. in this paper, we are concerned with young 
Cepheid-hosting clusters only, so the TP-AGB/AGB cor¬ 
rections are not critical. The Padova isochrones are com¬ 
puted over a large age range (from 6.60 to 10.10 dex in 
log(t/years)) which should comfortably capture all FM 
Cepheids. The Padova models use a moderate amount 


We fit the cluster CMDs using MATCH (Dolphin 


20021. We run MATCH in a mode that fits a single stel¬ 
lar population (SSP) to the cluster CMD. MATCH gen¬ 
erates model CMDs which include the effects of cluster- 
specific artificial star tests, and thus has some defense 
against spurious photometry effects in these crowded 
fields. The methodology, fitting parameters, and data 
products (photometry, artificial star tests, etc.) are the 
same as those used in the (forthcoming) CMD analy¬ 
sis of the full PHAT cluster sample (L. Beerman et al. 
in preparation). MATCH accounts for potential field 
star contamination in the cluster CMD by incorporat¬ 
ing a smooth binning of likely field stars situated in 
an annulus 10 times the size of the cluster aperture in 
the fit. We fix the cluster distance modulus at 24.47, 
and allow the metallicity of the cluster to vary from 
logio(^min) = -0.20 to logio(^max) = 0.10 in steps 
of 0.1 dex. The runs are performed with a resolution 
in logio(</years) of 0.10 over the range (6.6,10.1); and 
a resolution in Ay of 0.05 over the range (0.00,3.00). 
These ranges encompass all physically-reasonable solu¬ 
tions, assuming the Cepheids are cluster members. By 
fitting all models in this grid, and marginalizing over Ay 
and metallicity, we obtain a posterior probability density 
function (PDF) for the cluster age assuming the model 
we’ve adopted. 

We perform the fits with and without the Cepheid in¬ 
cluded in the data. This serves several purposes. First, 
isochrone fits are most robust (and degeneracy is mini¬ 
mized) when they are determined by main sequence stars 
rather than a handful of bright evolved stars. This is 
especially worth considering in these cases, where the 
Cepheid is often the only star in the corresponding sec¬ 
tion of the CMD. Secondly, this enables us to perform a 
meaningful test of the period-age relation. If the CMD 
fits are controlled by the Cepheid photometry, then the 
resulting MATCH age estimate is effectively a measure¬ 
ment of the Padova isochrone age at a given blue loop 
magnitude. Finally, since we are not certain that the 
Cepheid belongs to a dominant SSP as we have assumed, 
we can look for large changes in the PDF when the 
Cepheid is removed, which might indicate that this as¬ 
sumption fails. The reported ages used in Section ^ are 
those derived from the CMDs without the CepheicL^ 

In most cases, the PDFs with and without the Cepheid 
are very similar. Six clusters display differences in the 
best-fit age at our 0.10 dex resolution, and three of these 
show differences in the median age estimates. We can 
produce a rough approximation of the standard error 
of these estimates from half the 84-16 percentile differ¬ 
ence. Only one cluster (CC2) shows a difference in both 
best-fit and median age which exceeds the joint errors. 
The actual shift in both age estimates is only 0.20 dex. 
Several other PDFs show hints of a bimodal shift when 
the Cepheid is removed. Qualitatively, the lack of large 
shifts in the fitting results indicate that the MATCH fits 
are indeed capturing an underlying SSP which is roughly 
consistent with the Cepheids. 

The likelihood of the best-fit SSP model relative to 
the null model (no SSP, i.e. based only on the empir¬ 
ical field star CMD) is another test of fit quality. We 
re-run MATCH on each cluster, but constrain it to use 
only the smoothed background model derived from the 
























































PHAT/PS Cluster Cepheids 


7 


cluster annulus. We are concerned with the quantity 
A/ = /best - /background where / Is the ht value MATCH 
returns. This / is a Poisson-statistics analogue to 
and thus roughly corresponds to —21n(P); so the quan¬ 
tity ^ PbestZ-Pbackground informs us how much 

more likely the best fit found is relative to the null back¬ 
ground fit. This likelihood ratio ranges from 10^ to 10^® 
for our cluster fits without the Cepheid. The outliers at 
the low end of this quantity have correspondingly broad 
normalized PDFs, by necessity. We conclude that there 
are no catastrophic failures in the fitting results, and pro¬ 
ceed with the full cluster age PDFs (derived without the 
Cepheids) in hand. 


4. RESULTS 

We now use our period measurements and cluster age 
estimates to test the period-age relation for Cepheids in 
M31. This test requires a robust analysis of our MATCH 
results and a set of model period-age relations to compare 
with. 


4.1. Fitting Our Period-Age Data 

We adopt a probabilistic approach to fitting our data. 
This enables us to properly incorporate the full (non- 
Gaussian) probability distributions we computed using 
MATCH for the cluster ages, and to account for errors 
in the period determination from Kodric et _al.| (20131 
(substantial in two cases, as shown in fi'igureElb 

We fit a linear model to the period and age data: 
logio(age/years) = TOlogiQ(period/days)-|-&, with Gaus¬ 
sian scatter a in magnitude about the line. Rather than 
fit with m directly, we use 0 = arctan(TO) so that a fiat 
prior doesn’t favor steeper slopes. Given our data (rep¬ 
resented by D) and our SSP model for the clusters (H), 
we would like to determine the probability function for a 
linear slope and intercept: V{{9,b,a} \ D,n). Baye’s 
rule informs us that 

V{{0,b,a} I D,U)^r{D I {e,b,a},Il)Vi{e,b,a}) 

where V{{0,b,a}) represents our prior expectation for 
the fit parameters. We adopt a very generous prior (fiat 
over (—7r/2 < 0 < 0) and all possible intercepts). To 
evaluate the likelihood function ViD \ { 0, 6, cr} , H), we 
integrate the product of 1) the MATCH-provided PDFs, 
2) the Gaussian period PDFs, and 3) the modeled Gaus¬ 
sian variance about the line (0, b, cr) over a grid of points 
for each Gepheid. The product of these likelihoods for 
all the Cepheids yields the full likelihood function. We 
apply Markov chain Monte Carlo (MCMC) (using the 
package emcee: Foreman-Mackey et al. 2013) to sample 
the resultant posterior distribution. 

The MCMC fit is displayed in Figure ^ Using the 
16th/50th/84th percentiles of the MCMC chain, we 
obtain estimates for the slope and intercept of m = 
—0.69/g;2o, b = 8.38/g;|g; and for the scatter, a = 
0.25lg'2o. The constraints are broad, as can be seen in 
the wide range of the sampled lines in Figure This is 
a product of the broad age PDFs, which in turn indicate 
the difficulty of assigning an age to low-mass clusters at 
the distance of M31. 


4.2. Comparison and Discussion 


We compare our results to several period-age relations 
in the literature. It is important to note that the choice of 
evolutionary model used (in fitting clusters or modeling 
Cepheids) will in general affect the period-age relation, 
as it affects the age-lum inosity relation for stars near the 
instability strip. From Bono et al. (2005), we compare 
with the classical, fundamental-mode, relation, for Z = 
0.02 (solar metallit y). This relation was der ived using 
stellar models from Pietrinferni et al. (2004) which do 
not include convective core-overshooting. 

It is more difficult to find period-age relations derived 
using models which incorporate overshooting, such as the 
Padova models we used for CMD fitting. Thus, we de¬ 
rive our own semi-empirical period-age relation from the 
Padova models. Specifically, we use the MATCH util¬ 
ity fake to sample the Padova isochrones in the range 
6.6 < log(age/years) < 8.5, at a resolution of 0.05 dex. 
We adopt a sufficiently high SFR to populate all instabil¬ 
ity strip crossings. All other parameters (except initial 
distance modulus) are kep t identical to their state in the 
CMD fitting (Section 3.3). We select the models lying 


within the T pff, luminosity boundaries derived by Bono 


et al. (2005) (solar metallicity, non-canonical). We then 


utilize the theoretical p ulsation rela t ions f or fundamen¬ 
tal Cepheids derived by |Bono et al.| (|2000|) to relate the 
model temperatures, luminosities, ahd masses to pulsa- 
tional periods. Then a simple linear fit is performed to 
the resulting points in period-age space, providing the 
desired theoretical prediction. 

For further comparison, we locate two additio nal re¬ 
lations from the literature. Magnier et al. (1997) calcu¬ 
late a period-age relation in a semi-em pirical way: they 
use tracks from Schaller et al. (1992) which allow for 
overshooting, and an empirical pe riod-luminosity rela¬ 
tion found for M31 Cepheids. And Efremov (2003) fit a 
period-age relation empirically to a sample ot Cepheids 
from the LMC bar, uti lizing integrated col or results 
based upon models from Bertelli et al. (1994) which do 
incorporate moderate convective core overshooting; we 
quote their most probable relation. We plot these period- 
age relations in Figure alongside our MCMC fits, and 
compare the fit parameters in Table The median esti¬ 
mates agree reasonably well with the ensemble of other 
values given the large uncertainties in our fit. 

It is interesting to note that the largest disagreement 
in slope value is between our empirical fit and the Padova 
model predictions we derived, despite the fact that both 
rely on the same set of stellar models. To check whether 
this discrepancy is due to our choic e of pulsationa l mod- 
els, we also applied the approach of Magnier et al. (1997) 
and used empirical period-luminosity relations to relate 
the theoretical model luminosities to pulsational peri¬ 
ods and re-derive a Padova period-age prediction. The 
MATCH fake utility provides ma gnitudes in common 
filter systems for the model stars (Girardi et al.||2008|). 
We re-derive the period-age relation using several differ¬ 


ent period-lu minosity relations: two from Wagner-Kaiser 


et al. (2015) for the N IR HST filter FlfiOW and for the 


1-band; and two from Fouque et al. (2007) derived us¬ 
ing galactic G epheids in the V and 1 band s. The slopes 
derived using Wagner-Kaiser et al. (2015) are close to 
the shallow theoretical Padova results described above: 
m ~ —0.55 for both F160W and I. In contrast, the slopes 
found using Fouque et al. (2007) are closer to our steeper 
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Fig. 6. — The period-age relation fit for our cluster Cepheids. The 
blue violin plots represent the age PDFs produced with MATCH; 
their width at each age value represents the likelihood of that age 
from the MATCH SSP model. The period error is not marked on 
this plot (see Figure!^. The gray lines are random samples from 
the MCMC chain of nts. The line marked by ye llow triangles is 
the Padova period-age relation derived in Section |4.2| the red cir¬ 
cle lirie^e 2 resents_ym solar-metallicity period-age relation derived 
by|Bono et al.| l|2005|l; the orange-upper an d pink-lower tria ngles 
represent the period -aee relations found by |Efremov| (|2003|l and 
[Magnier et al.|||1997[l, respectively. 


TABLE 2 

Period-Age Relations 
logio(age/years) = mlogjQ(Period/days) -|- b 


Source 


m 

b 

Scatter 

M31 CC fit 
Padova fit 

o.Dy_Q 20 
-0.53 

Q OO+0.37 
8.38_q 30 

8.40 

y-^o_o.20 

0.07 

IBono et al.lfl 

2ooq 

-0.67 

8.31 

0.08 

iMaemer et al. 

(11997 ll 

-0.6 

8.4 

- 

|Etremov| 

Drol— ^ 

-0.65 

8.50 

- 


^ Predicted scatter in magnitudes abont the best-fit relation. 


empirical fit: m ~ —0.78, — 0.69 for V, I respective l y. We 
expect that the results from Wagner-Kaiser et al. (20151 
will be more robust, since their sample size is larger and 
Cepheid variability is reduced at long wavelengths. This 
reinforces the suggestion from the theoretical Padova fit 
that we should expect a shallower slope. Confirming this 
experimentally would require further testing on clusters 
with better constraints than our M31 sample can pro¬ 
vide. 

We conclude that our Cepheid age measurements are 
in broad agreement with existing theoretical and empir¬ 
ical period-age relations. There are several important 
caveats to note. Comparing ages derived using stellar 
models with differing assumptions about overshooting is 
not strictly valid, as overshooting parameters have an ef¬ 
fect on the derived age at fixed luminosity. In addition. 


we have compared period-age relations derived at differ¬ 
ent metallicities (solar and LMC/SMC). Our constraints 
on the period-age relation are too loose to distinguish 
these second-order effects; but they must be recognized 
when applying these tools to derive age estimates. 

5. CONCLUSION 

In this paper, we presented and analyzed a set of 
11 Cepheids in M31 which appear to be members of 
partially-resolved clusters. The sample was fo und by 
cross- c orrela ting the M31 Cepheid sample from |Kodric| 
et al. (12013) with the T HAT cluster sample presented 
by Johnson et af. (2015). We located the Cepheids in 
the HS'l' photometry by their variability and CMD lo¬ 
cations. This revealed that the PSl mean magnitudes 
were systematically brightened by crowding. We then 
performed fits to the selected cluster CMDs without the 
Cepheids. Our derived period-age relation for these ob¬ 
jects is poorly-constrained, but agrees broadly with the 
body of existing relations in the literature. 

The sample of cluster Cepheids we present are partic¬ 
ularly well-characterized for objects at the distance of 
M31. Combining the HST and PSl photometry allows 
the former to be properly dephased and the latter to be 
crowding-corrected. Our magnitude calibration results 
indicate that the PSl mean magnitudes are biased bright 
at the 1 magnitude level in the high est stellar density re¬ 
gions. Wagner-Kaiser et al. (|2015 ) uses PHAT photom¬ 
etry to further explore the eitect of such contamination 
on the period-luminosity relation for M31 Cepheids. 

Cluster Cepheids are essential calibrators for period- 
age relations and stellar evolution generally. This work 
constitutes the first direct test of the period-age rela¬ 
tion by CMD-fitting beyond the Galaxy and Magellanic 
clouds. Our results support application of such period- 
age relations to M31 and other solar-metallicity galaxies. 
Further investigation of the full PSl Cepheid sample in 
this context, once extended to more o f the M31 halo and 
to larger periods ( Kodric et al.||2013 ), could provide ad¬ 
ditional insight into the recent star formation history of 
M31. 
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Fig. 7. — Cepheid identification for CCl; see Figure[^for explanation of plot symbols. 



Fig. 8. — Cepheid identification for CCS; see Figure]^ for explanation of plot symbols. 


APPENDIX 

Here we present the diagnostic plots produced during Cepheid identification and dephasing described in Section 
Refer to the caption of Figurel^for a description of the various panels and symbols. In each case, the putative Cepheid 
is identified quantitatively by the variability metrics (bottom panels); and its HST measurements are compared to the 
PSl data (top panel). 
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Fig. 9. 
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Cepheid identification for CC4; see Figure]^ for explanation of plot symbols. 
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Fig. 10. — Cepheid identification for CCS; see Figure|^for explanation of plot symbols. 
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Fig. 11. — Cepheid identification for CC6; see Figure!^ for explanation of plot symbols. Note that this cluster was only observed 
PHAT field, and thus no variability information is available for the HST CMDs. Accordingly, the point size in the bottom panels is 
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Fig. 12. — Cepheid identification for CC7; see Figure|^for explanation of plot symbols. 
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Fig. 13. — Cepheid identification for CCS; see Figure|^for explanation of plot symbols. 
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Fig. 14. — Cepheid identification for CC9; see Figure|^for explanation of plot symbols. 
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Fig. 15. — Cepheid identification for CCIO; see Figure[^for explanation of plot symbols. 





Fig. 16. — Cepheid identification for CCll; see Figure[^for explanation of plot symbols. 
























